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After the PAMELA finding of an increasing positron fraction above 10 GeV, the 
experimental evidence of the presence of a new electron and positron spectral 
component in the cosmic ray zoo has been recently confirmed by Fermi-LAT. 
We show as a simple phenomenological model which assumes the presence 
of an electron and positron extra component peaked at ~ 1 TeV allows a 
consistent description of all available data sets. We then describe the most 
relevant astrophysical uncertainties which still prevent to determine source 
properties from those data and the perspectives of forthcoming experiments. 
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1. Introdution 

Recent experimental results raised a wide interest about the origin and the 
propagation of the leptonic component of the cosmic radiation. 

Among the most striking of those results, there is the observation per- 
formed by the PAMELA satellite experiment that the positron to electron 
fraction e+/(e^ + e+) rises with energy from 10 up to 100 GeV at least 
(Adriani et al. 2008^). This appeared in contrast with the predictions of 
the standard cosmic ray scenario and could therefore be interpreted as the 
smoking gun of new physics, unless a very soft electron spectrum was as- 
sumed. 

The significance of this anomaly increased when the Fermi-LAT space 
observatory measured the + spectrum in the 7 GeV - 1 TeV energy 
range with unprecedented accuracy and found it to be compatible with a 
power- law with index 7(e±) = -3.08 ±0.05 (Abdo et al. 2009,^ Ackermann 
et al. 2010"^); this slope is significantly harder than what estimated on the 
basis of previous measurements: the hypothesis of a steep spectrum was 
therefore excluded. 
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More recently, the same collaboration provided a ftirtlicr. and stronger, 
evidence of the positron anomaly by providing direct measurement of the 
absolute e+ and e~ spectra, and of their fraction, between 20 and 200 GeV 
using the Earth magnetic field. A steady rising of the positron fraction was 
observed by this experiment up to that energy in agreement with that found 
by PAMELA. In the same energy range, the e~ spectrum was fitted with 
a power-law with index 7(e~) = —3.19 ± 0.07 which is in agreement with 
what recently measured by PAMELA between 1 and 625 GeV (Adriani et 
al. 2011^). Most importantly, Fermi-LAT measured, for the first time, the 
e+ spectrum in the 20 - 200 GeV energy interval and showed it is fitted by 
a power-law with index 7(e+) = —2.77 ± 0.14. 

We will show in the following paragraph how all those measurements 
rule out the standard scenario in which the bulk of electrons reaching the 
Earth in the GeV - TeV energy range are originated by Supernova Rem- 
nants (SNRs) and only a small fraction of secondary positrons and elec- 
trons comes from the interaction of CR nuclei with the interstellar medium 
(ISM). Then we will see how the alternative scenario in which the presence 
of electron + positron component peaked at ~ 1 TeV is invoked allows 
a consistent description of all the available data sets. Finally we will dis- 
cuss to which extent astrophysical and particle physics uncertainties still 
aff'ect our modeling of cosmic ray leptons origin and propagation and how 
forthcoming measurements are expected to reduce those uncertainties. 

2. The necessity of a primary extra-component 

After the release of Fermi-LAT + e+ spectrum, it was clearly pointed 
out in several papers (see e.g. Grasso et al. 2009^ and Di Bernardo et al. 
2011^) that both Fermi-LAT and PAMELA measurements described in the 
Introduction are in contrast with a standard single-component scenario in 
which positrons are the secondary products of the nuclear component of 
cosmic rays (CRs) interacting with the interstellar medium (ISM). 

The main problems encountered by this kind of models can be summa- 
rized as follows. 

• As explained many times (see e.g. Serpico 2011*" for a recent re- 
view), they cannot reproduce the rising positron-to-electron ratio 
measured by PAMELA and recently confirmed by Fermi-LAT; 

• They are unable to reproduce all the features revealed by Fermi- 
LAT in the CRE spectrum, in particular the flattening observed 
at around 20 GeV and the softening at ~ 500 GeV. In fact, if 
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Fig. 1. Fermi-LAT and PAMELA data on electrons + positrons and electrons are com- 
pared to a double component phenomenological model. The absolute positron spectrum 
is compared to a single and double component phenomenological model. Red dotted 
line: e+ in single- component scenario. Red dot-dashed line: e+ in double- component 
scenario. Blue triple dotted- dashed line, black solid line: e~ and e~ + e+ in double- 
component scenario. Blue dashed line: e" diffuse background in double- component sce- 
nario. The Kolmogorov diffusion setup is adopted. 



such models are normalized against data in the 20 - 100 GeV en- 
ergy range, where systematical and theoretical uncertainties are the 
smallest, they clearly fail to match CRE Fermi-LAT and PAMELA 
e~ data outside that range. A different normalization results in even 
worse fits. 

With the release of the e~ and e"*" separate spectra by the Fermi-LAT 
collaboration the problems with the single component scenario became even 
worse. In fact, the e+ spectrum (Fig. 1) is clearly inconsistent with the pre- 
dictions of a single component scenario computed with DRAGON numerical 
diffusion package (and similar results are obtained with GALPROP). Even 
without considering numerical models, the simple consideration that the 
reported positron spectral slope is —2.77 ± 0.14 reveals how these data are 
incompatible with a purely secondary origin from proton spallation on inter- 
stellar gas: the source slope should be the same as the proton spectrum, i.e. 
~ —2.75 (Adriani et al. 2011®) and no room is then left for the unavoidable 
steepening due to energy- dependent diffusion and energy losses. 

A double component scenario is the most straightforward solution to 
these problems. 
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Fig. 2. Fermi-LAT and PAMELA data on the positron ratio are compared to a sin- 
gle and double component phenomenological model. Dot-dashed line: positron ratio in 
single- component scenario. Dotted line: positron ratio in double- component scenario 
due to conventional secondary positron production. Solid line: positron ratio in double- 
component scenario including extra- component. The progagation setup and modulation 
potential are the same of Fig. 1. The solar modulation potential is taken <I> = 550 MV 
in all figures of this paper. 

The idea dates back to the pioneering work by F. Aharonian and A. 
Atoyan 1995^ and was extensively studied after the release of ATIC and 
PAMELA data in 2008 (see e.g. Hooper et al. 2009^" and Profumo 2008"). 

More recently, we contributed to several papers in which it was shown 
that a consistent interpretation of the e"'"+e~ spectrum measured by Fermi- 
LAT and the PAMELA positron fraction can be naturally obtained in that 
framework (Grasso et al. 2009,^ Ackermann et al. 2010,'^ Di Bernardo et al. 
2011^). 

For example, in Fig. 1 and Fig. 2 we show that the double component 
model proposed in Ackermann et al. 2010'^ reproduces the data mentioned 
above and also the e"*" and e~ separate spectra, and their ratio, recently 
released by the Fermi-LAT collaboration and not yet available at the time. 
The model represented in those figures assumes a propagation setup charac- 
terized by a cylindrical diffusive halo with half-thikness of 4 kpc; a diffusion 
coefficient scaling with rigidity like p^^^ (corresponding to a Kolmogorov- 
like diffusion within the quasi-linear approximation) and a relatively strong 
reacceleration (the Alfven velocity is va = 30 kms^^). Solar modulation 
is treated here as charge independent in the force field approximation by 
fixing the modulation potential against proton data taken in the same 
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solar phase. In that model, the standard e~ primary component is tuned to 
fit Fermi-LAT data at low energy in the presence of the extra-component 
becoming dominant at higher energies; the injection slope for the primary 
electron component is set to —2.70 above 2 GeV, while under that energy 
a slope of —1.6 is adopted, in accord with recent constraints from the syn- 
chrotron spectra (see Jaffe et al. 2011^^). The extra component, instead, 
originates from a primary source of electron-t-positron pairs; it has an in- 
jection spectrum modelled in a simple way as a power-law with index —1.5 
plus an exponential cutoff at 1.2 TeV; the spatial distribution of this source 
is the same as the standard one and the propagation parameters are also the 
same; the normalization is tuned so that Fermi-LAT and PAMELA data 
at high energy are matched by the sum of standard + extra component. 
Both components are computed with DRAGON (even if it was checked that 
the same result can be obtained with GALPRDP). 

An issue remains open about the origin of the discrepancy between the 
prediction of this, or similar, models and the positron fraction measured 
by PAMELA below 10 GeV. In the next section we will show as that dis- 
crepancy may be interpreted as the consequence of an incorrect choice of 
the propagation setup and discuss other uncertainties which can affect the 
electron and positron spectra in that low energy range. 

3. LOW ENERGY. Impact of astrophysical uncertainties 




Fig. 3. Effect of changing the diffusion halo height. Solid line: h = 1 kpc; dashed: h = 
10 kpc. 
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Fig. 4. Effect of the diffusion setup. Solid line: KRA; dashed line: KOL. 



Cosmic ray electrons and positrons, either belonging to the standard or 
the extra component, propagate in the Galaxy undergoing several physical 
processes: diffusion, reacceleration, energy losses. Such complex motion is 
effectively described by a well-known diffusion-loss equation (Berezinskii 
et al. 1990^'^). In this equation several free parameters are involved: the 
height of the halo in which the propagation takes place, the normalization 
and energy dependence of the diffusion coefficient (the latter parametrized 
by the parameter 5), the Alfven velocity that influences the effectiveness 
of reacceleration; moreover, several astrophysical inputs need to be consid- 
ered: the injection spectrum, the spatial distribution of the source term, 
the interstellar radiation field, the gas distribution. 

The free parameters that appear in the diffusion-loss equation are 
constrained by some CR observables such as Boron-to-Carbon (B/C) or 
antiproton-to-proton ratio; different diffusion setups exist in the literature, 
obtained through comparison of experimental data with the prediction of 
semi-analytical codes (Maurin et al. 2001,^^ Donato et al. 2004^^) or numer- 
ical packages such as DRAGON or GALPROP (see e.g. Di Bernardo et al. 2010^^ 
for DRAGON- related models and Trotta et al. 2011^'' for a GALPRDP-based 
analysis) . 

The uncertainties related to the diffusion model and to the astrophysical 
inputs were discussed in the latest years in several papers making use of 
semi-analytic codes (e.g. Delahaye et al. 2010^®). In the following we will 
briefly analyse the impact of these uncertainties adopting the DRAGON code. 
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One of the most relevant parameter is the halo height. According to the 
analytical computations by Bulanov and Dogel/^ while at low energy the 
electrons (or positrons) are distributed throughout all the diffusion halo, as 
the energy increases the electrons occupy a smaller and smaller fraction of 
the halo due to energy losses. This is relevant especially for the secondary 
positron spectrum. In fact, since their injection power is determined by the 
CR nuclei density in the Galactic disk, a thicker halo results in a larger 
dilution of their density in the halo hence a in smaller flux on the Earth. 
Numerical computations confirm the expectation of this heuristic argumcmt 
as shown in Fig. 3. From the plot it is also evident that large halo heights 
are disfavoured by the data. 

Even fixing the height of the diffusion halo, the choice of the diffusion 
setup can also affect the low energy spectra of CR leptons. This is evident 
from Fig. 4 where we compare the predictions of two different models which 
both reproduce nuclear CR data: 

• a Kraichnan-like diffusion setup with ^ = 0.5 and moderate reac- 

cclcration (that was pointed out as the preferred one in a DRAGON- 
based maximum likelihood analysis with focus on both B/C and 
antiproton high energy data^^) 

• a Kolmogorov-like diffusion setup with 5 = 0.33 and high reaccelcr- 
ation (that was pointed out as the preferred one in a GALPROP-based 
maximum likelihood analysis with focus on B/C data^^) 

It is clear from that plot that the Kraichnan-like setups allows a better 
fit of low-energy positron ratio measured by PAMELA; this consideration, 
together with several other facts (high reacceleration models do not permit 
a good fit of antiproton data and cannot reproduce the spectrum of the 
synchrotron emission of the Galaxy), led us to conclude that models with 
strong reacceleration are disfavoured. 

4. High energy uncertainties and the nature of the 
extra-component 

In the double component scenario discussed in Sec. 2, the positron spectrum 
above ~ 10 GeV is dominated by the primary extra component. The nature 
of its source is one of the hottest matter of debate in the CR physics. 

Galactic pulsars were suggested as natural source candidates of a pri- 
mary CR positron component well before PAMELA results (Aharonian and 
Atoyan, 1995.^) More recently, it was noticed that a single, nearby, pulsar 
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(such as Monogem or Gcminga) could explain the positrons fraction excess 
found by PAMELA (Hooper et al. 2009^°). 

In the Fermi-LAT era, we showed (Grasso et al. 2009^ and Di Bernardo 
et al. 2010*^) that also the 6+ +e~ measured by that experiment can consis- 
tently be explained in the same terms: if one considers the observed nearby 
pulsars within 2 kpc and assumes that a relevant fraction of their rota- 
tional energy is transferred into e+ + e~ pairs (~ 30%), under reasonable 
assumpions on the injection spectrum and cutoff it is possible to reproduce 
all existing data. In the cosmic ray channel, this scenario has two possible 
testable consequences: 

• the detection of a CR electron anisotropy towards the most relevant 
sources (in our analysis, Monogem and Geminga^'^); 

• the presence of some bumpiness in the e~ and e"*" spectra in the 
TeV region due to the contribution of several pulsars. 

Those two signatures are somehow complementary: if a single pulsar 
give the dominant contribution to the extra component a large anisotropy 
and a small bumpiness should be expected; if several pulsars contribute the 
opposite scenario is expected. 

So far no positive detection of CRE anisotropy was reported by the 
Fermi-LAT collaboration, but some stringent upper limits were published. 
In Di Bernardo et al. 2010^ we showed that the pulsar scenario is still 
compatible with these upper limits. Also, no evidence of spectral bumpiness 
has been found so far in the e"*" -|- e~ spectrum. 

It should be noted that several astrophysical uncertainties prevent ac- 
curate predictions of the CRE anisotropy and of the spectral bumpiness. 
For example, unknown irregularities in the local structure of the Galactic 
magnetic field may distort the angular distribution of the CRE flux due to a 
nearby pulsar. Furthermore, due to the stochastic nature of the emission 
of nearby SNRs, the CRE standard component is expected to be subject to 
fluctuations which may produce anisotropics and spectral bumpiness which 
may hide those due to pulsars. 

The other possible scenario to explain the origin of the extra component 
is more exotic but very appealing as it invokes DM annihilihation/decay as 
the origin of the extra component. Plenty of papers were published on 
that subject after the release of PAMELA and Fermi-LAT results (see e.g. 
He 2009^° for a review). That scenario, however, present some problems. 
The most important are the following ones. 

• It requires a heavy DM particle mass - O(TeV) - and an annihi- 
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lation cross section much higher than that predicted by standard 
cosmology if one assumes that DM is a thermal relic. 
• Since no excess was detected for antiprotons, the annihila- 
tion/decay channels must include only leptons (lepto-philic DM). 

Although several DM models which may fulfil those conditions were 
developed, another issue arises when electroweak corrections are taken into 
account. Those corrections, in fact, give rise - even in a Icpto-philic scenario 
- to soft electroweak gauge bosons, and hence to antiprotons, at the end 
of their decay chains (Ciafaloni et al. 2011^^). Since those exotic p are 
produced mainly in the Galactic Center region, the flux reaching the Earth 
strongly depends on the properties of CR propagation in the Galaxy. As we 
discussed in Sec. 3, these properties are still subject to strong uncertainties. 
It was shown in Evoli et al. 2011^^ that, accounting for those uncertainties, 
a scenario in which a heavy DM particle annihilates into muons is still 
compatible with the antiproton constraints. In the same paper it was also 
shown that AMS-02 is expected to constrain even more these models since 
its sensitivity to antiprotons will be much higher. 



5. Conclusions and future perspectives 

In this contribution we argued as recent experimental data rule out the 
standard scenario in which CR positrons are produced only by CR spal- 
lation onto the ISM and showed as an empirical model which invokes an 
extra component fulfils all data sets. We also discussed several uncer- 
tainties which still prevent to infer some of the properties of CR electron 
and positron sources. We argued that at low energy those uncertainties 
are dominated by our poor knowledge of CR propagation (which prevent 
an accurate determination of the injection spectrum of the e~ standard 
component) while at high energy the effect of the stochastic nature of as- 
trophysical sources prevails (which makes more difficult to decide between 
the astrophysical and DM origin of the extra component). 

Forthcoming experiments like AMS-02 and CALET are cixpcctcd to rc;- 
duce drastically the uncertainties on the propagation parameters by provid- 
ing more accurate measurements of the spectra of the nuclear components of 
CR. Fermi-LAT and those experiments are also expected to provide more 
accurate measurements of the CRE spectrum and anisotropy looking for 
features which may give a clue of the nature of the extra component. 
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